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Cross-Kerr nonlinearity in optomechanical systems 
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We consider the response of a nanomechanical resonator interacting with an electromagnetic cavity 
via a radiation pressure coupling and a cross-Kerr coupling. Using a mean field approach we solve 
the dynamics of the system, and show the different corrections coming from the radiation pressure 
and the cross-Kerr effect to the usually considered linearized dynamics. 


I. INTRODUCTION 

Cavity optomechanics offers a framework to study the 
coupling between an electromagnetic field and the vibra¬ 
tions of a mechanical resonator. The interaction between 
these two systems is usually mediated by a radiation- 
pressure type coupling proportional - through a coupling 
constant g— to the number of photons n c in the cavity and 
the displacement of the mechanical resonator. The radi¬ 
ation pressure coupling offers the possibility of altering 
the resonant frequency of the mechanical resonator and 
its damping. The latter can be used for cooling @H5] or 
amplification [7j. Moreover, the nonlinearity of the in¬ 
teraction may allow for the observation of macroscopic 
quantum phenomena such as quantum superpositon of 
states mm or quantum squeezed states J5|. The re¬ 
quirements for observing these quantum phenomena are 
the necessity of being close to the ground state and be¬ 
ing in the strong coupling regime mm where g is larger 
than the cavity and the mechanical resonator decay rate. 
However, g is usually weak and to bypass this constraint 
a strong drive to the cavity is applied at the cost of losing 
the nonlinear property of the interaction. 

Our recent proposal HD], in which the cavity and the 
resonator are coupled to a Josephson junction, shows that 
the interaction between the cavity and the resonator can 
be enhanced via the non-linearity of the Josephson effect. 
Quadratic and higher-order interactions in the displace¬ 
ment have been investigated also in different setups such 
as membrane in the middle geometries [HHT3]. 

Analogously to the setup mentioned above, the non¬ 
linearity of the Josephson effect leads to an additional 
nonlinear interaction, namely a cross-Kerr coupling g c k 
between the cavity and the resonator. The difference 
resides in the fact that in the Josephson junction setup 
the relative value of g c k and g depends on the value of 
the gate charge to a superconducting island, whereas in 
mm, it generally reflects the position of the resonator 
within the cavity. 

In the context of optomechanical systems, the cross- 
Kerr coupling between the resonator and the cavity in¬ 


duces a change to the refractive index of the cavity de¬ 
pending on the number of phonons in the resonator, 
whereas the radiation pressure coupling gives rise to an 
analogous effect, but depending on the displacement of 
the mechanical resonator. 

In this paper we solve the dynamics of the cavity and 
the mechanical resonator in the presence of the cross- 
Kerr and the radiation pressure coupling. We determine 
the effects of the cross-Kerr coupling on the red and blue 
sidebands within a mean field approach. In particular, 
we demonstrate that the sideband peak is shifted due 
to the cross-Kerr coupling. In addition, the cross-Kerr 
coupling induces a nonmonotonuous response of the ef¬ 
fective mechanical damping as a function of the number 
of photons pumped into the cavity. 
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FIG. 1. Schematic picture of the system. A cavity and a 
mechanical resonator coupled via a radiation type coupling 
g and a cross-Kerr coupling g c k . The number of photons in 
the cavity n c is coupled to the oscillations of the mechanical 
resonator x and the number of phonons in the mechanical 
resonator n m - 


II. MEAN FIELD APPROACH 
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We consider an electromagnetic cavity with frequency 
w c and linewidth k coupled to a mechanical resonator 
with frequency and linewidth 7 . The number of 
phonons in the cavity n c is coupled to the vibration am- 
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plitude of the mechanical resonator x via a radiation- 
pressure type coupling g. In addition the number of pho¬ 
tons n c is coupled to the number of phonons n m in the 
mechanical resonator through a cross-Kerr coupling g ck 
(Fig-0- The Hamiltonian of the system is (H = 1) 

H = uj c a' a + ujmb^b — ga l a( b^ + b) — gckoJab'b, (1) 

where a and b are the annihilation operators of the cavity 
and the mechanical resonator, respectively. We treat the 
interactions with a mean-field (MF) approach. Within 
this frame, the radiation-pressure interaction becomes 


ga^a^ +b) = g ((a^)a + (a)a^ — (a^a))(b^ + b) 

+ (a) a — (a)a^ — (a*)a)(b^ + b) , ( 2 ) 


where (A) stands for the average of A over the static 
nonequilibrium state of the system (mean field). The 
negative terms in Eq. Q are included to suppress dou¬ 
ble counting. The first line of Eq. ([2| describes exchange 
processes between the resonator and the cavity while the 
second line gives a frequency shift of the cavity which 
is proportional to the average displacement of the res¬ 
onator. This decomposition allows us to find the usual 
results of the weak radiation pressure coupling [7) [16] . In 
MF, the cross-Kerr coupling becomes 


g ck a)atfb = g ck 


(o t a)6 t 6+ (6 t 6)a t a 


+ ( 6 t a) 6 a t + ( 6 a t ) 6 t a + (ba)tfa ] + ( 6 t a t ) 6 a 


■ (3) 


The term (a^a) 6^6 ((&!&} era) describes a Hartree-like in¬ 
teraction between the resonator and the cavity while the 
other terms describe exchange processes between the res¬ 
onator and the cavity. Thus we can rewrite the Hamilto¬ 
nian as 


H =[u c - g ck (b ] b)]a)a + [u m - g ck (a^a)}tfb 
— G[(a^)ab^ + (a)aW] — G*[(a^)ab + (a)a' i b] 

+ g[(a) a — (a)a^ — (a , )a)(b^ + b) — (a^a)( 6 ^ + 6 )], (4) 

where the expectation values of the different operators 
have to be determined self-consistently within the MF 
picture and G = g + g c k(b)- We assume the usual experi¬ 
mental situation where w c and where the cavity is 

driven with a coherent field of strength f p oscillating at 
frequency oj p = lo c + A. Using the input-output formal¬ 
ism mi the equations of motion are 

a = -*[— A - g ck (tfb)\a - j k/ p 

+ iG*(a)b + iG(a)b 1 — ig(b 1 + b)[a — (a)] (5) 

b = -i[uj m - g c k(a < a)]b - -b + ^ 7 bin 

+iG(a l )a + iG(a)a* — ig(a l a), ( 6 ) 


Here we have written the cavity operator a in a frame 
rotating with frequency ui p neglecting the fast oscillati- 
ing terms. We define bi n to be the thermal input of 
the resonator satisfying (fr» n (i)) = 0 and ( b\ n (t)bi n (t')) = 
n th 5(t — t'), where n th is the number of phonons in the 
thermal bath damping the resonator. We split the cav¬ 
ity and the mechanical operators into a sum of coherent 
and fluctuation parts, i.e., a = 6a + a and b = 6b + /3 
with a = (a), /3 = ( b ) and (6a) = (6b) = 0. As usual, 
we assume that a and /? oscillate at the same frequency 
as the coherent drive so that a = $ = 0. With these 
approximations, the solutions of Eqs. § (§) are 

f 

“ = f-i[A- 5cfc (6t6)-(G*/3 + G/3*)]’ '' 1 

= i(2G-g)\a\ 2 -ig(8a^5a) 

'y/2 + i(u} m -g c k(a^a))) 


In the derivation of Eqs. |7| | 8 ), we have assumed, in 
agreement with what is usually done in the optomechan¬ 
ical literature (see e.g. [IS]), A + g(tf + b) « A. The 
equations of motion for the fluctuations in the Fourier 
space are given by 


~K ./ X \" 

~-'(uj + A) 
2 - i(u - Wm) 


6a 

5b 


iGaStf + iG*a5b (9) 

iGa* 5a + iGaSa I + y/^b in (10) 


where A = A + gck^b) and C6 m = co m — g ck (a i a). The 
effect of the thermal drive 6 j n on the response of the 
cavity is mediated by the coupling G. Through this cou¬ 
pling the oscillations of the mechanical resonator produce 
sideband peaks at u>d ± in the cavity response. They 
allow for the exchange of energy between the cavity and 
the resonator when A w ±ui m mi HZ]- These processes 
are depicted in Fig. [2] For A m —u: m the system is in the 
red sideband regime and one can transfer energy from the 
resonator to the cavity, thus the mechanical resonator is 
damped and cooled. For A « u> m , the system is in the 
blue sideband regime and one can transfer energy from 
the cavity to the resonator, thus the mechanical resonator 
is excited and heated. In order to find the correction to 
the damping, we solve the response function of 8a for 
the thermal input 8bi n . We find that it is a Lorentzian 
function peaked at w m + u> s hift with 


(^shift — 


|G| 2 |a| 2 (A 2 -^+^) 




+ (“m + A) 2 V + (d 



and whose linewidth is 7 + r opt with 


( 11 ) 


r opt = \g\M z k 


1 


“t + (6j m + A)* 



( 12 ) 
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FIG. 2. Cooling (heating) process. The cavity is driven with 
a frequency uid. = to c — cn m {oJd. = w c + w m ). The drive does 
not allow a transition from |n m , n c ) —¥ \n m , n c + 1) but allow 
the transition from \n m ,n c ) — > |n m — 1 ,n c + 1) (|n m ,n c ) —>• 
| n m + 1 ,n c + 1)). The cavity relaxes then to the state \n m — 
l,n c ) (| n m + 1 ,n c )) resulting into cooling (heating) of the 
mechanical resonator. 


Integrating the Lorentzian function obtained above, we 
obtain the number of phonons and photons coming from 
the thermal vibrations of the resonator. We get [ 6 ] 


(Sb^sb) = inth + T °P tnm ° , 

'y “i - r 0 pt 

(5a) 5a) = G 2 \a\ 2 (Stf5b) 

y j + (w m + A)" 


(13) 


(14) 


1 


^ + A) 2 


with 


— 


(A m + A ) 2 + 
4A u) m 


(15) 


Eqs. Q, Q, (13) and (14) form a set of self-consistency 
equations and are solved in the next sections in order to 
find the number of phonons in the resonator and photons 
in the cavity. We now focus on the sidebands. 


III. OPTIMAL COOLING/HEATING 


In order to minimize/maximize the optical damping 
r op t, we set A = T&m ■ The upper sign refers to the 
red sideband (T opt > 0 ) and the lower sign to the blue 
sideband (T opt < 0). In the resolved sideband limit, 
ojm !?> k 7 , the frequency shift and the optical damp¬ 
ing become 


^shift — “F 


T op t = ± 


|G | 5 


\a\ 


= =F- 


|G | 5 


a 


Wm bJm dcki^®) 

| 2|„,|2 


4|G| 2 M : 


K 


(16) 


(17) 


The result for the optical damping Eq. ( |T7| ), is identical 
to the one usually obtained in optomechanics in the ab¬ 
sence of the cross-Kerr coupling. The effect of g c k shows 



^shift bk >0 

FIG. 3. Schematic picture of the red sideband with and with¬ 
out cross-Kerr coupling for A < 0. For g c k > 0 the sideband 
peak is shifted to lower values while for g c k < 0 the sideband 
peak is shifted to higher values. 


up only in the frequency shift which now depends on the 
number of coherent and thermal photons in the cavity 
Eq. (16). In Fig. [3] we show a schematic picture of what 
happens to the sideband in the presence of the cross-Kerr 
coupling. 

In the Doppler limit (w m < n) the frequency shift and 
optical damping are given by 


Wshift = T4|G| 2 |a| 2 ^ 

T 

r _i_ 4|G | 2 |ck | 2 


^m Qck (G Cl) 

+ 4 (uj m ~ g c k(a f a)) 2 
4(uj m - g c k(a f a)) 2 
+ 4 (u) m - g c k(a f a)) 2 


(18) 

(19) 


Now both the frequency shift and the optical damping 
depend on the cross-Kerr coupling. In Figs. [4][5]we plot 
the number of phonons and the optical damping as a 
function of w m /« for the red sideband in the Doppler 
limit. Since the cross-Kerr coupling shifts the mechanical 
frequency, the value of r opt is shifted as well. The sign of 
the shift is given by the sign of g c k■ Otherwise we recover 
the cooling of the resonator for the red sideband (Fig. [4]) 
and the parametric instability when r opt = —7 for the 
blue sideband (Fig. [ 5 ]). 


IV. CASE WITH A=uj m . 


In experiments the parameter one can tune directly is 
the detuning A and not A as it can be difficult to set 
A = u} m for each value of a as the pump strength is 
varied. Therefore, another regime we consider is the case 
where A = yu m , i.e, setting A = uj m + g c k(b^b). In 
this case the frequency shift and optical damping in the 
red (upper sign) and in the blue (lower sign) sideband 
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FIG. 4. Steady-state phonon number in the resonator and 
F op t (inset) as a function of the ratio uJ m /n at the red sideband 
for the optimal case A = — Cj m with 7 = 10” 3 k, g = 1 CP 2 k. 
The number of photons pumped into the cavity is fixed to 
| O' | 2 = 100 and the bath temperature corresponds to n th = 
100 . 



FIG. 5. Steady-state phonon number in the resonator and 
F 0 pt (inset) as a function of the ratio a; m /fv at the blue 
sideband for the optimal case A = Cjm with 7 = 10 _ 2 k, 
g = 10 - 2 k. The number of photons pumped into the cavity 
is fixed to \a\ 2 = 100 and the bath temperature corresponds 
to n th = 10. The dashed lines indicate the onset of the para¬ 
metric instability for r opt = — 7 . 



FIG. 6 . Steady-state phonon number in the resonator and the 
optical damping F op t (inset) as a function of the number of 
photons pumped into the cavity in the case where A = —w m . 
The values for the parameters are 7 = 10~ 3 u; m , k = 10 - 1 ai m , 
g = 10 - 2 w m and the bath temperature corresponds to n th = 
100 . 


become 

^shift 


r opt 

In Figs. [6] and [7] the steady-state phonon number and 
the optical damping are plotted as a function of the 
number of photons pumped into the cavity, for the red 
and blue sidebands respectively. For the red sideband 
(Fig. | 6 j) the optical damping increases with increasing a. 
When g 2 k \a\ 2 k 2 /4 the optical damping becomes in¬ 

versely proportional to the number of photons pumped 
into the cavity, consequently, the cooling deteriorates 
when pumping more phonons in the cavity. 

In the blue sideband (Fig. [ 7 ]) the main effect of a small 
cross-Kerr coupling is to limit the instability to a finite 
number of phonons, (tfb) ~ \f^h\G\\a\/\g ck \ + |a| 2 . 
This thus competes with the usual limitation coming 
from the intrinsic (Duffing) nonlinearity of the resonator. 

In Figs.(8|9]we plot the frequency shift as a function of 
the number of photons pumped into the cavity for the red 
and blue sidebands. For the red sideband when g ck > 0 
( g ck < 0 ) the frequency shift increases (decreases) as a 
increases until a ~ (m b) after which it increases (de¬ 
creases). For the blue sideband when g ck > 0 the fre¬ 
quency shift decreases while for for g ck < 0 it increases. 
The difference at small a between the red and blue side¬ 
bands arises from the fact that in the red sideband when 
pumping more photons into the cavity the cooling im¬ 
proves, thus the number of phonons in the mechanical 


t IG| 2 H 2 

[k 2 /4 — 2 g ck u m {{tfb) - (a^q))] , , 

[g 2 ck (m-(a^a)r + K y 4 ] ’ 1 ; 

= ± _ l G l 2 H 2 * _ (21) 

9c k ((^b) — (ata)) 2 + k 2 /4’ 1 J 
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a) 



b) 



FIG. 7. a) Steady-state phonon number in the resonator and 
b) the optical damping r opt as a function of the number of 
photons pumped into the cavity in the case where A = oj m . 
The values for the parameters are 7 = 10 _ 3 u; m , k = 10 _ 1 w m , 
g — 10~ 2 uOm and the bath temperature corresponds to n th = 
10 . 


resonator decreases, making it possible to have a num¬ 
ber of photons in the cavity of the same order and larger 
than the number of phonons in the resonator. 


V. CONCLUSION 

In conclusion, we have solved the dynamics of a me¬ 
chanical resonator coupled to an electromagnetic cavity 
via a radiation pressure coupling and a cross-Kerr cou¬ 
pling using a mean field approach. We have shown that 
the cross-Kerr coupling shifts the frequency of the me¬ 
chanical resonator and of the optical cavity, the shift 


depending on the number of photons in the cavity and 
phonons in the resonator. In addition, we have shown 
that when the detuning of the pump is equal to the fre¬ 
quency of the mechanical resonator the variation of the 
optical damping saturates instead of being linearly de¬ 
pendent on the number of phonons pumped into the cav¬ 
ity. 



FIG. 8 . Frequency shift as a function of the number of pho¬ 
tons pumped into the cavity when A = — oj m with 7 = 
10 ~ 3 LUm k = 10 ^ 1 oi m , g = 10~ 2 uOm and the bath tempera¬ 
ture corresponds to n th = 100 . 



FIG. 9. Frequency shift as a function of the number of pho¬ 
tons pumped into the cavity when A = ca m with 7 = 10 _ 3 a) m 
k = 10 _ 1 o; m , g = 10~ 2 oj m and the bath temperature corre¬ 
sponds to n th = 100 . 
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